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We show that the known matrix representations of the stationary state algebra of the 
Asymmetric Simple Exclusion Process (ASEP) can be interpreted combinatorially as vari- 
ous weighted lattice paths. This interpretation enables us to use the constant term method 
(CTM) and bijective combinatorial methods to express many forms of the ASEP normali- 
sation factor in terms of Ballot numbers. One particular lattice path representation shows 
that the coefficients in the recurrence relation for the ASEP correlation functions are also 
Ballot numbers. 

Additionally, the CTM has a strong combinatorial connection which leads to a new 
"canonical" lattice path representation and to the "w-expansion" which provides a uniform 
approach to computing the asymptotic behaviour in the various phases of the ASEP. 

The path representations enable the ASEP normalisation factor to be seen as the partition 
function of a more general polymer chain model having a two parameter interaction with a 
surface. 
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1 Background and notation. 



The Asymmetric Simple Exclusion Process (ASEP) is a simple hard core hopping particle model. 
It consists of a line segment with r sites. Particles are allowed to hop on to site 1 if it is empty, 
with rate a. Any particles on sites 1 to r — 1 hop on to a site to their right if it is empty, with 
rate 1. A particle on site r hops off with rate /3 - as illustrated in figure ^ The state of the 
system at any time is defined by the set of indicator variables {ri, • • • , Tr}, where 

1 if site i is occupied 
otherwise 

and the probability, P{t; s) of the system being in state r = (ri, • • • , r^), at time s given some 
initial state at s = 0, satisfies a master equation (see j4j for details). 
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Figure 1: The ASEP model. 

For a review of the extensive literature on this model see Derrida pp. A significant step 
forward in the understanding of the mathematical aspects of the model was made with the 
realization, by Derrida et al 01, that a stationary state solution, Ps{t), of the master equation 
could be determined by the matrix product Ansatz 

r 

Ps{r) = {W\ HinD + (1 - n)E)\V)/Z2r (1.1) 

i=l 

with normalisation factor given by 

Z2r = {W\{DEY\V) (1.2) 
provided that D and E satisfy the DEHP algebra 

D + E = DE (1.3a) 

where {W\ and \V) are the eigenvectors 

{W\E = -{W\, D\V) = hv). (1.3b) 
a p 



2 



These equations are sufficient to determine Psi'^) but Derrida et al [1] also gave several inter- 
esting matrix representations of D and E and the vectors \V) and {W\, any one of which may 
be used to determine -Ps'(t). 

The primary result of this paper is to show that each of the three matrix representations of 
the DEHP algebra can be interpreted as a transfer matrix for a different weighted lattice path 
problem. This then allows the normalisation, correlation function and other properties of the 
ASEP model to be interpreted combinatorially as certain weighted lattice path configuration 
sums - see section 2. One of the path connections is similar to that discussed in Derrida et. al. 

The lattice path interpretation has two primary consequences: The first is that it provides a 
starting point for a new method (the "constant term" or CT method) for calculating the normal- 
isation and correlation functions - see section 3. This reproduces several existing results (but 
by a new method) and also provides several new results. One of note, the "w-expansion" , arises 
from a rearrangement of the constant term expression which leads to form of the normalisation 
in terms of the variables uJc = a(l — a) and tOd = /3(1 — /?). The coefficients in this expansion are 
Catalan numbers, the asymptotic form of which enables a uniform approach to computing the 
asymptotic behaviour of Z2r as r ^ oo in the various phases of the ASEP model - see section 
5. The results agree with those found in j^, by steepest descent methods. 

A bonus following from the lattice path interpretation of the algebra representations is that 
one of the lattice path interpretations (a slight variation of representation 3) has a natural 
interpretation as a polymer chain having a two parameter {ni,K2) interaction with a surface. 
In this context is a partition function for the "two-contact" polymer model - see section 
4. We also obtain recurrence relations (on the length variable) for the partition function of this 
polymer model and hence also for the ASEP normalisation. 

The second primary consequence of the lattice path interpretation follows from the CT 
method itself, as the CT method has very natural combinatorial interpretations. For example, 
the normalisation can be written in several different polynomial forms depending on which 
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variables you use: 

Z2r = Y,P^n}na''l3"' (see - g^l) (1.4) 

n,m 

Z2r = Y.Pnl'^"d"^ (see- 633) (1.5) 

n,m 

Z2r = Y.PnL^>2 (see - g^Hl) (1-6) 

n,m 

Z2r = E^^Sl^^^"2" (see - g^U) (1-7) 

n,m 

where all the polynomial coefficients, pn,m, are integers. (All the variables in these polynomials 
are related by simple equations eg. a = 1/a, c = a — 1 etc. - see above and in later sections). 
Since, in each of these cases, the normalisation arises from a weighted lattice configuration sum, 
all the above coefficients have a direct combinatorial interpretation as enumerating a particular 
subset of the paths (eg. those with exactly m steps with the first weight and exactly n steps 
with the second weight). 

However, we show that each of the above polynomial coefficients has an alternative combi- 
natorial interpretation which corresponds to enumerating a different, unweighted, set of lattice 
paths eg. in the coefficient, Pm]n was given as 

(1) m(2r — m — 1)! 

Pn,m-n = r\{r - m)] 

which, with a simple rearrangement, can be seen to be a particular "Ballot number" 

- R 

where Bg^h enumerates Ballot paths (see section 3.1 for details) of length s and height h. Thus 
Pn^rn-n which enumerates a special set of paths with m weights of type a and n — m weights of 
type /? = is seen to be determined in terms of the much simpler combinatorial problem of 
enumerating unweighted Ballot paths of length 2r — m — 1 and height m—1. This correspondence 
between between the two combinatorial problems (one pair for each coefficient) arises as a 
bijection between the two path problems. This result may be turned around: If a bijection 
between a particular Ballot path problem and the weighted path problem can be proved then it 
provides an alternative derivation of the normalisation polynomial. 

Finally, in section 6, the recurrence relations for the ASEP correlation functions derived in 
[1] are shown to follow from the lattice path interpretations. The coefficients of the terms in the 
recurrence relation are also seen to be various Ballot numbers. 
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2 Matrix Representations and Lattice Path Transfer Matrices. 



Derrida et al provided three different matrix representations for tlie ASEP algebra. Repre- 
sentation one, 



wliere 



representation two, 



fp P P P P ■■■\ 

1 1 1 1 •• 

1 1 1 •• 

1 1 •• 

\l \ \ \ I 

(M^il = (1, a, c?, 6? , . . . 

a = 1/a, 



/i 1 --A 

1 1 •• 

1 1 •• 

1 1 •• 

Vi ; ; ; j 



{W2\ = k{1,c,c',c^...) 

where 

c = a — 1. 

and representation three 

fp K --A 



110 

110 

1 1 

y; ; ; ; 



^1 



fo --A 

1 •• 

1 •• 

1 •• 

v; ; ; ; 



|yi) = (i,o,o,o,...f 

P = i/P, 

/l 
110 
E2= 110 
110 

\\ \ I \ 

\V2) = K{l,d,cf,d^...)' 



d = P-l, 
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(2.1) 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



(VFsl = (1,0,0,0,...) 11/3) = (1,0,0,0,...)^ 
5 



(2.6) 



where 



= a + P — aP = 1 — cd. 

Since these matrices and vectors satisfy the algebraic relations ()1.3() , the normalisation factor 
H1.2() for the ASEP model, can be evaluated using any of the three formulae 



'2r 



{Wj\{DjEjY\Vj), J = 1,2,3. 



(2.7) 




D 



E 



2 4 6 

3 
5 
7 



13 5 7 
/ A \ 

2 
4 
6 



Figure 2: a) Action of D and E matrices on the square lattice, h) the labelling of the matrix 
elements. 

Each of these three representations can be interpreted as the transfer matrix of a particular 
weighted lattice path problem. If the rows of the Dj matrix and the columns of the Ej matrix 
are labelled with odd integers 'Lo^d = {1)3,5,...} and the columns of Dj and rows of Ej are 
labelled with even integers l^even = {0, 2, 4, ...}, then {Dj)^/ is the weight of a step from an odd 
height k to even height i and [Ej)k^i is the weight of a step from an even height j to an odd 
height Since the rows and columns are labelled with non-negative integers the steps are only 
in the upper half of 1? . 

Similarly, the elements of {Wj\ and \Vj) are labelled by Zf^^^ and are the weights attached to 
the initial and final vertices of the paths. The matrices Dj and Ej act successively to the left 
on the initial vector (yVj\ and {Wj\{DjEjY\Vj) is the weighted sum over all paths of length 2r 
which begin and end at odd height above the x— axis. 

An example path for each of the three representations is shown in figure 13 Notice that for 
the first representation the paths with non-zero weight begin at any odd height and end at unit 
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height, for the second they both begin and end at any odd height and for the third they begin 
and end at unit height. All these paths are defined explicitly below. 




Figure 3: Examples of the paths for each of the three representations a) representation one, 
"Jump step paths", h) representation two, "Cross paths" and c) representation three, "One up 
paths". 

Derrida et. al. gave another expression for the normalisation factor (see equation (39) of 
This form does not arise directly from any of the above three representations, however, we will 
show (corollary that it is the partition function , defined in (|2.1(ij) . corresponding to a 
"canonical" path representation (see figureElfor an example). In ^01 we provide a combinatorial 
derivation of the equivalence of the above three and a number of other path representations of 
the normalisation factor to the "canonical" representation. 

2.1 The Lattice Path Definitions. 

We consider paths whose steps are between the vertices of the half plane square lattice H = 
{(x,y)| X £ Z,y £ where Z (resp. Z+) is the set of integers (resp. non-negative integers). 

Definition 1 (Lattice paths). A lattice path, oj, of length t > is a sequence of vertices 
(vo, vi, . . . ,vt) where Vi = {xi,yi) G H, and for t > 0, Vi — G 5f where {5f , i = 1, . . . ,t} is 
the set of allowed steps. For a particular path, u), denote the corresponding sequence of steps by 
£■{00) = 6162 ■■■Ct, Cj = (fj_i,i;j). A subpath of length k of a lattice path, uj, is a path defined 
by a subsequence of adjacent vertices, (fj, Uj+i, . . . , Uj_|_fc), ofuj. 

Definition 2 (Dyck P^^\ Ballot p/.f\ & Cross ^^^2 Paths). The set of length 2r 
Dyck paths, P2^\' length t and height h Ballot paths , and length t and heights hi and /12 
Cross paths P^^^^ all have step set S,f = {(1, 1), (1, — 1)}. Dyck paths have vq = (0,0) and 
V2r = (0, 0), Ballot paths have vq = (0, 0) and vt = {t, h) and Cross paths have vq = (0, hi) and 
vt = {tM), hiM > 0. Denote, P^? = {Jn^Mez^,, P^r-MM 

7 



a) Dyck path 



b) Ballot path 



c) Cross path 




Figure 4: An example of a a) a Dyck path, b) a Ballot path, c) a Cross path, d) an Anchored 
Cross path, e) a Hovering path, f) a Separated Hovering path, g) a One Up path, h) a Jump 
Step path and i) a Marked separated hovering path. 

Definition 3 (Elevated Dyck path (Bubble)). An elevated Dyck path or Bubble, is a 

subpath, {vi, . . . , Vi^k), k > 0, for which yi = yiJ^k o-nd yj > y-i for alli<j<i + k, k>0. If 
k = then the elevated Dyck path is a single vertex. 

Definition 4 (Anchored Cross Pj:!^fl^\^ paths). The subset P^"^^\^ C Pj^^l of Anchored 

(c) 

Cross Paths is defined as all the paths in P^.f^^^ with at least one vertex in common with the 
line y = l. Denote, P^f^ = [Jh.M&^aa ^^MM 

Definition 5 (One Up p!^^ paths). The set of One Up paths, p!^^ is the set of lattice 

paths of length 2r which have step set Sf = {(1, 1), (1, —1)}, vq = (0, 1) and V2r = (2r, 1), i.e. 

p(0) _ p(C) 
^2r — -'^2r;l,l- 

Definition 6 (Jump Step P^*;^^ paths). The set of Jump Step paths of length t, Pj^.f^ , is the 
set of lattice paths which have step set 



5, 



J 



{(1, —1)} for i even 

{(1,-1)} U {(1,2^ + 1)|£ G Z+} fori odd 



with vq = {0, h), h G "^odd, "^t = {t, 1) . The "height", gi, of a step = (vj-i, fj) is defined as 
Hi ~ Vi-i- Odd steps with g-i £ {3, 5, . . . } will be called "jump" steps. If gi = 1 (resp. gi = —1 ) 
the step is called an "up" (resp. "down") step. Denote, ^2/^ = Uhez ^^^'Zr'^h 

Note: Jump step paths never visit the x-axis since such a visit can only occur on an odd 
step and return to y = 1 is impossible since all even steps are down. 

Definition 7 (Hovering PgV ■> Separated Hovering 2p Marked Separated Hov- 
ering P2^2p paths). Hovering paths, P2^^ are One Up paths with no vertex on the x-axis. 
Separated Hovering paths, P2r^l, < p < r, are the subset of the Hovering paths which have 
V2p = (2p, 1). The vertex V2p is marked (with an empty circle - see figure^ ) and known as the 
separating vertex. Marked separated hovering paths, P^2p ' ^'^^ obtained from P2r^l by 
marking (with an solid circle - see figure^ ) subsets of the steps (or vertices) which return to 
y = l. 

Note: For all paths considered Xj + yi is either odd for all i or even for all i, i.e. the paths 
are confined either to the odd sublattice or the even sublattice. 

Definition 8 (Contacts, Returns). A vertex of a Ballot or Dyck path in common with the x- 
axis is called a contact. All contacts except the initial one are called returns. The polynomial 
Rt{h;K) = ^ (s) kP^^\ where p{uj) is the number of returns for the path ui, is called the 
return polynomial for Ballot paths. 

2.2 Weights and lattice path representations. 

For each of the three different representations the following lemma converts the matrix formula 
1)2. 7() for the normalization factor into a sum over one of the path sets defined above, where the 
summand is a product of the step weights u)j*^^(ej), an initial vertex weight Wj{vo), and a final 
vertex weight vuj {vt). 

(i) 

Lemma 1. The normalisation factor, Z)^^ for each of the three matrix representations, j = 
1,2,3, can be written as 

Zg)= (2-8) 



with 

2r 



i=l 



wUv2r) (2.9) 



9 



where Pg^^ = P2r \ -^ir^ — ^2r'^ > '^'^^ — ^2^^ ' ^''^^ weight W^^\u)) of a particular path 
Lo with step sequence, £{co) = ei . . . e2r is defined, for each of the three cases, as follows. 
For j = 1 



for j 



w{{({),2k + l)) =a^ keZ+ 
' _ 

Pifei = {{i - 1, 1), (i, 2k)) , k e Z+ and i odd 
1 otherwise 

U)i{v2r) = 1, 

wi{{0,2k + l)) = Kc'' keZ+ 



w. 



step / 



W 



2 (ei) = 1 
^((0,2A: + 1)) =Kd'' 



k G 



and for j = S 



^i((o,i)) 



w 



step / 



K i/e, = ((i-l,l),(i,2)) or a = {{i- 1,2), {i,l)) 

p if a = {{i- 1,1), {i,o)) 

a ifei = {{i-l,0),{i,l)) 

1 otherwise 



wk{o,i)) = 



(2.10a) 
(2.10b) 

(2.10c) 



(2.11a) 
(2.11b) 
(2.11c) 



(2.12a) 



(2.12b) 



(2.12c) 



Proof. The above lemma is a direct consequence of the lattice path interpretation of the Dj and 



Ej matrices as transfer matrices. 



□ 



The equivalence of the different expressions is a consequence of the invariancc of the nor- 
mahsation factor under similarity transformations relating the different matrix representations 
of D and E. 

In order to compute Zg,, , it turns out to be a little more convenient to rearrange the weights 
associated with representation three. If we do so we obtain the following corollary. 
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Figure 5: An example of a representation three, one up path, with the weights re-organised 



(3) 

Corollary 1. An equivalent set of weights for Z^', is 



w;i((0,l)) = l 



w 



step I 



Ki = a(3 i/ci = ((z - 1,0), (i, 1)) 
= k2 if(>. = ((^_ l,2),(i,l)) 

1 otherwise 



^{((0,1)) = 1 



(2.13) 



(2.14) 



(2.15) 



(3) 

Note: this rearrangement of the weights is only valid for computing Z2,/ ■ For correlation 
functions one has less freedom in the weight rearrangement. 




2s 



Figure 6: An example of a canonical path - a Separated Hovering path with a VF^^^ weighting. 

Later in the paper (corollary we will show that the normalisation factor can also be 
expressed in terms of Separated Hovering paths. 



(2.16) 



p=0 



where W^^^ {ujp) is defined by with 



wimi)) = l 



-u;^*^P(ei) = < 



a if Ci = ((i - l,2),(i,l)) and i < 2p 
(3 ifei = ((i-l,2),(i,l)) and«>2p 
1 otherwise 



"^^((0,1)) = 1 



(2.17a) 



(2.17b) 



(2.17c) 
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Thus, for any particular path, Wp, ah the a. weighted steps (if any) occur to the left of vertex 
(2p, 1) and all the /5 weighted steps (if any) occur to the right of (2p, 1). We call this combina- 
tion of paths and weights the "canonical" path representation of the normalisation factor. An 
example is shown in figure El 

Lemma 2. Lei P'^.2y be the set of marked separated hovering paths obtained from P2r^p by 
marking subsets of the steps which return to y = 1 then 



^^? = E E ^^'"^(-.) (2-18) 

where the weight W^'^'^^{ujp) has a factor c for each marked return step which occurs to the the 
left of V2p and a factor d for the other marked return steps. 



The lemma is proved in ^0] in two stages. First an involution on -P2^^ for which p!>°'^^ 



2r 



IS 



the fixed set and then a bijection between ^2^*^^ Up=o ^2t-2p ■ 
Corollary 1. 

7(2) _ 7(5) 

a=l+c,/3=l+d 

Proof. Substituting a = l-|-c, /9 = l-|-(iin the weight attached to a path cjp G -^2*5 ™^ 
expanding leads to a sum of terms obtained by weighting each return step to the right of v^p 
(if any) with either 1 or d and the remaining returns (if any) by 1 or c. Marking the subsets of 
steps weighted c or d determines a set of paths belonging to P'^^p ■ □ 

3 Methods. 

In this section we briefly review methods to be used and results obtained previously [7] as they 
will be required in the next section. One additional new lemma is stated. 

3.1 The constant term method. 

Definition 9. The constant term operation, CT[-], is defined by 

CT[/(2;)] = the constant term in the Laurent expansion of f{z) about z = 0. 

The number of t— step lattice paths with step set S^^ which begin at (0, 0) and end at {t, y) 
with no further constraint (i.e. replacing the constraint y £ Z"^ in the definition of H by y G Z) 
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is the binomial coefficient {i(^^_yy) for wliich the constant term formula is 



CT[{z + l/zYzy]. (3.1) 



t 

By the reflection principle |Sj, the number of t— step Ballot paths of height h is obtained by 
subtracting the number of unrestricted paths which begin at (0, —2) from those beginning at 
(0, 0), both ending at (t, h). 

Bt,h = \pI^\ = CT[{z + l/zfz''] - CT[{z + l/zYz^+^] = CT[A'z^{l - z^)] (3.2) 

where A = z + 1/z. Differencing the binomial coefficients expresses Bf fi in terms of factorials; 
for t + h even 

{h + l)t\ 

{\{t + h) + m\{t-h))v ^'-'^ 

Dyck paths are Ballot paths ending at (t, 0) so the number Dyck paths with 2r steps is 



IP^V^^I = Cr[A*(l - z')] = J = a, (3.4) 
a Catalan number. 

Instead of using the reflection principle the Ballot numbers may be obtained as the solution 
of the equations, t,h > 1, 

Bo,h = Shfi (3.5a) 

Btfi = Bt-i,i (3.5b) 

Bt,h = Bt-i^h-i + Bt-i^h+i, (3.5c) 

all of which are solved by CT[A*z''(l - z^)]. 

Historically, Ballot numbers arise in the combinatorial problem of a two candidate election. 
If you ask how many ways can t votes be caste such that the first candidate ends h votes ahead 
of the second candidate and at any stage of the voting never has fewer votes that the second 
candidate. 

3.2 Partition function for the one contact model. 

Previously |7j we proved the following proposition concerning the return polynomial (see defini- 
tion |H1). 
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Proposition 1. The return polynomial for Ballot paths of length t and height h is given by 

ii,(/.;K) = cr[^-A__^] (3.6) 

Proof. For t,h > 1, Rt{h; k) is determined by the recurrence relations 

Ro{h; k) = 5h,o (3.7a) 
Rt{0;K) = KRt^i{l,K) (3.7b) 
Rt{h; k) = Rt-i{h - 1; k) + Rt-i{h + 1; k). (3.7c) 

The first and last equations are the same as for the unweighted Ballot paths and are satisfied 
by CT[A*2;'^(1 — z'^)g{z)] provided that on expansion g{z) has no negative powers and g{0) = 1 
(which will be the case). We have introduced the factor g{z) to allow the second equation to be 
satisfied. But z + 1/ z = kz -\- {1 — {k — l)z'^)/z and choosing g{z) = 1/(1 — (k — l)z'^) 

CT[K\l - z^)g{z)\ = CT[A*-i(z + l/z)(l-z2)^(^)] 

= KCT[k'-^z{l - z^)g{z)] + C7r[A*-i(l/z - z)]. 

The result follows since the last term may be evaluated to give zero. □ 

Corollary 2 ([7j). 

i(t-h) i(t-h) 

Rt{h-K)= Bt^h+2ra{'^-ir = Bt.m-l,m+h-l^"' ■ (3.8) 

m=0 m=0 

Proof. The first equality follows by expanding (|3.6() in powers of {k — l)z^ and using ()3.2() . The 
second is obtained by rewriting the constant term formula as 

Rt{h- k) = CT[ 4—^] (3.9) 

I — KZ/A 

and then expanding in powers of kz/A. This result suggests a bijection between Ballot paths 
of length t with m returns and Ballot paths of length t — m — 1 and height m + h — 1 and 
hence suggests a combinatorial proof. Such a bijection was given in ^.9 . A Ballot path can 
be represented schematically as shown in figure [7| The Bubbles (see definition represent a, 
possibly empty, arbitrary elevated Dyck path. Using this schematic representation, proving (|3.8|) 
is then straightforward. For simplicity we show the bijection in the case /i = in figure |H1 The 
expansion in powers of k — 1 was also obtained by bijection in 0. In this case the bijection is 
between Ballot paths with at least m returns, m of which are marked, and Ballot paths of the 
same length but of height h + 2m. □ 
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Figure 7: a) An example of a Ballot path showing the "terraces" which define the Bubbles shown 
schematically in b) which represents of a Ballot path of length t and height h = 4. 

a) . b) c) 




2j *■ ■* 2r— m »■ -< 2r— m— 1 



Figure 8: a) Schematic representation of a Dyck path with m = 3 returns, b) All return steps 
deleted c) First step deleted produces a Ballot path of height m — 1 and length 2r — m — 1. 

^) ^ m marked returns b) 




Figure 9: a) An example of a marked return Ballot path of height h = 1 with m = 2 marked 
returns and b) its schematic representation in terms of bubbles, c) Each marked return rotated 
through 90 degrees anti- clockwise (or equivalently replaced by an up step) producing a Ballot path 
of height 2m + h and length t. 

The above methodology is typical of that used for the more complicated two parameter case 
in the next section. A constant term formula will be derived and then rewritten in four different 
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ways each giving rise to an expansion in pairs of different variables the coefficients of which are 
Ballot numbers and are shown, by bijection, to enumerate various types of lattice path. 
In the next section we will also need the following lemma. 



Lemma 3. Let P2r^l be the set of Separated Hovering paths (definition^) th 



en 



^ l-^2r;2il — l-^2r+2l (3.10) 



s=0 

Proof. Since we can shift any Hovering path down on to the x-axis to give a Dyck path with 
vertex (2s, 0) marked, the lemma says the number of Dyck paths with one contact marked is 
equal to the number of (unmarked) Dyck paths two steps longer. A sketch of the bijective 
combinatorial proof is shown in figure 1101 The construction is a bijection, since, given any 
length 2r + 2 Dyck path a unique length 2r marked contact Dyck path is determined by deleting 
the rightmost step and the rightmost step from y = to y = 1 (and marking the left vertex of 
the latter step). □ 

a) b) c) 





o r o 

marked contact inserted edge at mark added edge 

Figure 10: a) Dyck (or hovering) path with one marked contact, b) insert an up step at the mark 
produces a unique height one Ballot path (length 2r+l) and c) adding a final down step produces 
a unique Dyck path of length 2r + 2 



4 The two contact model. 

We begin by computing using the W^'^'^ weights - we will refer to this as the two contact 
model. We will generalise the model by allowing arbitrary starting and ending heights for the 
paths. In particular, let the paths be of length t, start at (0,y*) with y' G '^odd and terminate 
at (t,y^), t + y-f e Zodd, i-e. Cross paths in P^^]^yf ■ 

This model can also be thought of as a polymer model with two different surface interactions, 
or "contacts" . The path interacts with the "thickened" surface y = 0, 1 via two parameters, ni 
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and K2 and the partition function is defined by 

Zt{yf\y';Ki,K2)= (4-1) 



t:y^,yJ 



where the sum is over Cross paths of length t with given initial and final heights and y^^ and 
the weig ht VF(^)(cj) is defined by (|2.9|) and (|2.13l) . This is a generalisation of the ASEP partition 



function Z^'^ thus 



Z^l^ = Z2r{l\l;aP,K^) (4.2) 



This generalisation allows this partition function to be determined by recurrence relations 
similar to those for the return polynomial of Ballot paths. By considering paths of length t — 1 
which can reach the point {t,y) by adding one more step the partition function Zt{y-^\y'^] ki, K2) 
may be seen to satisfy the equations, 

Zi(l|2/^Ki,^j2) = (4.3) 

Zoiy\y';Ki,K2) = 5yy (4.4) 

Zt{0\y';Ki,f,2) = Zt-iil\y';Ki,K2) (4.5) 

Zt{l\y';Ki,K2) = KiZt-iiO\y';Ki,K2) + K2Zt-ii2\y';Ki,K2), t>2 (4.6) 

and for t = 1, 2, . . . ; y = 2, 3, . . . and t + y odd 

Zt{y\y'] Ki, K2) = Zt-i{y - l\y'; ki, K2) + Zt-i{y + ki, ^2)- (4.7) 

These partial difference equations can be solved for Zt{y^y^ ; ki , K2) by using the constant 
term method. 

Proposition 2. With Ri = Ki — \, z = 1/ z and K = z + z, for y* G '^odd o-nd y > 1 

Zt{y\y'; ^1, ^2) = Cr[A*(z^-^' - z^+^"-2)] + (5,,,i + ^2(1 - 6y.^,))CT[zy+y'-^A'G{z)] (4.8) 
where 



1- {ki+ K2)Z^ - K2Z^ 



Note: For y* = and y > 1 

Zt{y\0;Ki,K2) = KiZt-i{y\l]Ki,K2) (4.10) 
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Proof. Substituting (|4.8|) into the partial difference equations and noting that G{z) may be 
expanded in even powers of z with no inverse powers shows that (|4.3|) . ()4.4() and ()4.7() are satisfied. 
H4.5() may be taken as the definition of Zt{0\y^; ki, K2) and H4.6() may then be transformed into 

Zt(l|y*; Ki, K2) = KiZt-2il\y';Ki,K2) + K2Zt-i{2\y';Ki,K2), t > 2 (4.11) 

To verify that this equation is also satisfied we note that 

A'^G{z) = A{l/z -z) + (ki + K2zA)Giz) (4.12) 

In the case = 1 this gives 

Zt{l\l; Ki,K2) = CT[A'G{z)] = CT[A'~\l/z - z)] + KiZt^2{l\y'; ^^l,^^2) + K2Zt^^{2\y'; ki,K2) 



and (|4.1H) follows since the first term is zero. 
Otherwise = 3, 5, . . . in which case 

Z(_i(2|y^ Aci, AC2) = GT[{z''-y' - + K2CT[zy' A'-^Giz)] 

and using this together with (|4.12() and the fact that the first term in (|4.8() vanishes when y = 1 

Zt{l\y';Ki,K2) = K2GT[zy'-^A'G{z)] 

= K2GT[zy'-\l/z - z)A*-i] + KiZt-2W; Ki, K2) + KlGT[zy'A'-^Giz)] 

= K2CT[izy'-^ - il/z)y'-^)A'-^] + K,Zt^2W; ^1, '^2) + ti2Zt^i{2\y'; ^1,^2) 

and again ()4.11|) follows since the first term evaluates to zero. 

An alternative constructive proof of this proposition is given in appendix B. □ 

Corollary 3. The number of t— step Cross paths which begin at = (0,/ii) and end at = 
(t, /12) is given by 

z,(.*;i.i) = (.,,_,'^^, J -(,,,_,;_,_,_ J (4.13) 

in terms of which we can write, Zt{y\y^; ki, H2) cls 

Ztiy\y'; k^, K2) = Zt{y - 2\y' - 2; 1, 1) + {5y.^^ + K2{1 - <5^,,i))Zf (y|y^ aci, K2) (4.14) 

where 

Z^iyW- ^1,^2) = CT [A'zy+y'-^G{z)\ (4.15) 
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is the partition function restricted to Anchored Cross paths except that for > 1 the step leading 
to the first visit to y = 1 has weight 1. 



In particular 

Ztiy\l;Ki,K2) = Z^(y\l;KuK2)=CT [A*G(z)z^-i] 



(4.16) 



Notes: 



It follows from the constant term formula that 

Z^{y\y';Ki,K2) = Zt{y + y' - 1|1;ki,K2) 



(4.17) 



• Setting Ki = K2 = 1 gives the number of unweighted paths starting at height 1. 

Zt{y\l; 1, 1) = CT[K\l - z^^-^) = CT\k'^\\ - z^)zy] = Bt+i,y 
as expected since the paths biject to Ballot paths by adding an initial up step. 
Proof. With Ki = K2 = 1, G{z) = 1 — and substituting in 1)4. 8() gives 

Zt{y\y^; 1, 1) = CT [A^zV-y' - zy+y'+') 
which yields (|4.13;) . 

Zt(t/ — 2|y* — 2; 1, 1) is the number of Cross paths which avoid y = 1. This follows since such 
paths are in simple bijection with the paths starting at height — 2 and ending at height y — 2, 



(4.18) 



eg. just push the whole path down (or up) two units. The second term in (|4.14|) is therefore the 
partition function for Anchored Cross paths. In the case > 1 Anchored Cross paths always 
have a first visit to y = 1 and the step leading to this visit has a factor /t2. Removing this factor 
leaves Z^{y\y^; ki, K2) which is therefore the partition function for Anchored Cross paths except 
that for y^ > 1 the step leading to the first visit to y = 1 has weight 1. Setting = 1 in (|4.14l) 
gives ()4.16|) since the first term vanishes. 



Corollary 4. 

Z^{y\y';KuK2) = CT 



(1 - z^)A^-^zy+y'-^ 

1 - {ki +zkK2){z/Kf 



□ 



(4.19) 



E 

i=o 



{t-y~y')/2-j+l 

E 



3 



k=0 



'j + A;^ 
, k J 



B 



t~2j-k-l,y+y'+k-3 



(4.20) 
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iVoie; When y = = 1 and j = t/2, k = the coefficient is indeterminate but setting 

_i = 1 gives the correct answer. This corresponds to the path which alternates between 
y = and y = 1. Notice that replacing the factorials in the definition of B^ h by Gamma 
functions gives Bt,t = 1 for all t > —1. 

Proof. Rewrite G{z) in the form 

1 - z2 

"^^^^ " (zA)(1-(ki + zAk2)(z/A)2)' 

expand in powers of ki and K2 and use the CT formula (|3.2|) for Ballot numbers. 

Equation (|4.2U|) may also be proved by a combinatorial argument. To simplify the proof we 
only consider the ASEP case y = = 1. The extension to general y, y* is straightforward. 
Recall that, from the weight definition, H2.13() . the returns to y = 1 are weighted with ki or H2 
depending on whether the return is from below or from above y = 1. The binomial coefficient 
corresponds to choosing a particular sequence of ki and k,2 weighted returns. For each particular 
sequence of returns we need to show there are Bt-2j-k-i,k-i possible path configurations. 

We first represent a particular sequence schematically and then show any path corresponding 
to the schematic can be bijected to a Ballot path with the correct height and length. Schemat- 
ically an example of a particular sequence of ki and K2 returns is shown in figure ITTI 




Figure 11: Schematic representation of a one-up path corresponding to the return sequence 

k\k2K,iK2- 

We now perform three operations to biject a given sequence into a Ballot path. 

• First delete all K2 return steps, see figure W^K ) . This produces a path of length t — k and 
height k + 1. 

• Next, delete the first up step above y = 1 (if any - which is the case if /c > 0), see figure 
112b ). This produces a path of length t — k — 1 and height k. 

• Finally, delete all 2j steps originally below y = 1, see figure ^J:). This produces a Ballot 
path of length t — k — 1 — 2j and height A: — 1 as required. 
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Given the sequence of ki and K2 the reverse direction for the bijection is obtained by simply 
reversing the forward mapping - see figure IT^ □ 




Figure 12: Bijection of first example sequence in figure [771 to a Ballot path: a) Delete all K2 
return steps, then h) delete first up step above y = 1, then c) delete all steps originally below 
y = 1, to give d) the final Ballot path. 

Corollary 5. 

(t-y-y^)/2+l {t-y-y')/2-j + l / ' + k\ 

Z^{y\y'-tii,K2)= Yl ''i'^M'^. ) Bt+k+i^y+y.+2j+3k-i (4.21) 

i=o fe=o ^ ^ 

Proof. Rewrite G{z) in the form 



1 — (ki + zAk2)z^ ' 

expand in powers of ki and R.2 and use the CT formula ()3.2() for Ballot numbers. 

As with corollary 01 the result, (|4.2H1 . may also be proved combinatorially as follows. Again, 
in order to simplify the proof we consider only the ASEP case y = y^ = 1. The substitution 
Ki = 1 + Ki means that a return which was weighted with k^ is now weighted with either Ri or 
1. The factor 

"'^ ^ Bt+k+i,2j+3k+i (4.22) 

is the number of paths with a subset of exactly j of the returns from below marked and exactly 
k of the returns from above marked corresponding to the choice of weight Ri or 1. The bino- 
mial coefficient in (|4.22j) is the number of ways of choosing a particular sequence of Ri and R2 
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weighted marks (reading the path from left to right), whilst, for a given sequence, the Ballot 
number represents the number of paths corresponding to the sequence. The most general path 

a) b) c) 




d) ^ e) O f) ^ 




Figure 13: a) Schematic One Up path or "frying pan" with no steps marked, b) Schematic One 
Up path with only the last step Ri marked, c) Schematic One Up path with only the last step 
R2 marked, d ) An example showing how the schematic frying pan represents a, possibly empty, 
One Up path, e) An example showing a frying pan followed by a down step and an up step 
(which forms a Ri marked return), f) An example showing a frying pan followed by an up step, 
then a Bubble then a final down step (which forms a R2 marked return). 

corresponding to a given sequence can be represented schematically by concatenating the cor- 
responding schematic sub-paths shown in figure [T3k ) and lLSb l with a final "frying pan" shown 
in figure Et). Examples of sub-paths corresponding to the three types of schematics are illus- 
trated in figures EJl) - I13t "). Note, the shaded regions of the schematics represent any number 
(possibly zero) of steps. Examples of two possible sequences are illustrated in figure im 

Thus, for a given return sequence we need to show that there are -Bt+fc+i.2j+3A:+i return 
marked paths. Without loss of generality we choose a typical sequence and represent it schemat- 
ically as shown in figure ITU 

Thus to prove the Ballot number factor in (|4.22|) we need to biject any schematic marked 
return sequence to a Ballot path of length t + k + 1 and height 2j + 3A; + 1. We do this by 
bijecting each schematic in the sequence to a sub-Ballot path (plus, possibly, an extra step) and 
then concatenate them all together. 

• Thus, the last frying pan, of length say, 2r', bijects to a height 1, length 2r' + 1 Ballot 
path - see figure El 
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a) 



K2 



1^2 



b) 




Figure 14: Schematic ki - K2 marked hovering path, a) For the sequence K2K1K1K2 and b) for 
the sequence RiRiR2Riiii- 



a) b) c) 




added edge 



Figure 15: Bijection of a "frying pan" to a height one Ballot path. The final path is one step 
longer than the original. 

• A marked K2 schematic of length 2ri bijects to a height 2, Ballot path with an additional 
final up step (see figure [T6|) . The final length is 2ri + 1 as an extra step has to be added. 

a) b) 




Figure 16: Bijection of a K2 schematic to a height two Ballot path, with a final up step. The 
final path is one step longer than the original. 

• Finally, a marked Ri schematic of length 2r2 bijects to a height 1, Ballot path with an 
additional final up step (see figure ITTjl . The final length is unchanged. 
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b) c) 

o 



move edge to front 

Figure 17: Bijection of a Ri schematic to a height one Ballot path, with a final up step. 

Putting these moves all together is illustrated in figure IT51 which show clearly a Ballot path of 
length t + k + 1 and height 2j + 2k -\- k + 1 is obtained. 

□ 

An explicit evaluation of the partition function of the canonical path representation, defined 
by (|2.16|) . is the case = = 1 of the following result. 

Corollary 6. With ki = a(3 and K2 = ol + (3 — a/? 

Ztiy'\y^;Ki,K2) = Yl 5t-m-l,m+y>+3;/-3 (4-23) 

m=0 " ^ 



m 



^(t-y'-yf)+l 

m=0 j=0 



Notes: 



• For the ASEP model {y^ = y^ = l,t = 2r) with r > 1, (|4.24|1 reduces to the result of 
Derrida et al equation (39) 

Z„.(1|1;m,k.) = t ""r^V' E^^"-' = (4.25) 
^-^ r'.ir — m)\ ^-^ 

m=l ^ ' j=0 

which they obtained directly from the algebra (|1.3|1 (0] appendix Al, equation (A12)) 
without the use of a matrix representation. 

• The result for the ASEP model may be written in terms of the return polynomial for 
Ballot paths, thus 

r r-j 



Z2.(1|1;Ki,K2) = X]«'Z]^2r-,"^-l-?'+^-l^' 
j=0 £=0 



J2R2r-jij-JW. (4.26) 

j=0 
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a) 



b) 



move edge move edge 




t - 
o 



c) 












^ added edge 









,i±2i. 



added edge 




t + k + 1 



Figure 18: Bijection between a One Up path with a given sequence of marked Ri and R2 steps 
and a Ballot path, a) The marked One Up path, b) After the application of move illustrated in 
figure \T1[ c) After the application of move illustrated in figure UR e) Final Ballot path after 
rotating the remaining marks. 

This formula may also be derived using a path representation based on recurrence relations 
of Derrida, Domany and Mukamel PI- 

Proof. By definition of a and (3 



1 - (ki + R2)z'^ - K2Z^ = [1 - (a - l)z\l -{(3- l)z2] = z2(A - az){A - pz) 



(4.27) 
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and from H4.9jl 



G(.) = ^. ''f ^ = I ^—^ I (4.28) 



{A-az){A- (3z) a -(3 \l-az/A l-fiz/k^ 

The result follows by expanding in powers of a and /3, substituting in ()4.15|) and using the CT 
formula (|3.2() for Ballot numbers. 

In the ASEP case = = l,t = 2r the coefficient -B2r-m-i,m-i in (|4.24|) is equal to the 
number of Dyck paths with m returns to y = (see (|3.8j) ). The equality of Z2r(l|l; ^1,^2) with 
defined by (ITTHl) . follows by raising the Dyck paths so that they become hovering paths 
u G P-z^"^ with returns to y = 1. The factor a^(5^~^ in ()4.24|) corresponds to weighting the 
first j returns of oj with a and the remainder with (3. The sum over s in p. 16)1 is obtained 
by partitioning the weighted paths according to the position (2s, 1) of the j*^ return (ie. the 
separation vertex). 

The equality of z'^^ and Zjj!^ is also shown directly in [IF by involution. □ 
Corollary 7. Let a = c + 1 and (3 = d + 1 then ki = (c + l)(d + 1), K2 = 1 — cd and 

{t-y-y*+2)/2 m 

m=0 i=0 
{t-y-t+2)/2 ^^^^^ _ ^^^^^ 

= Y Bt^l,y+yi+2m~l ^ (4-30) 

771=0 

in particular, in the ASEP case 

r m 

Z2.(l|l; (c + l){d + !),!- cd) = Zf^ = B2r+i,2m+iYcid"'~^ (4.31) 

m=0 j=0 

Note: Equation 1)4. 31() is not given in 11] but (34) and (35) of I3J together give the related 
formula 

I — i J — i 

This expression involves an infinite series whereas our expression is finite. By corollary |21 the 

coefficient in ()4.32() is the number of Cross paths of length 2r with hi = 2i — 1 and /12 = 2j — 1 

(C) 

and the double sum extends over • '^^^ factor 1 — cd restricts the sum to Anchored Cross 
paths. The equivalence of ()4.32|) and 1)4. 31|) is shown in ^Uj by constructing an involution on 
Pg^^ having Pg"'^^'' is its fixed point set. 
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Proof. From (|07|) 

1 - (ki + R.2)z 
and substitution in (|4.9I) gives 

G{Z): 



K2Z 



(l-cz2)(l-dz^ 



(1 - z^)zA 



(l-Cz2)(l-(iz2) 

The result follows from (|4.15|) using the expansion 

1 1 / c 



d 



(4.33) 



(4.34) 



{1 - cz^){l - dz^) c-d\l-cz^ l-(iz2 
expanding in powers of z and using the CT formula (|3.2|) for Ballot numbers.. 

Again a combinatorial proof is possible which for simplicity we only give in the ASEP case. 

(2) 

The equality with Zj^, follows from corollary^ To obtain the Ballot number formula we obtain a 
bijection between (Js=o ^2r-2s^ °^ Ballot paths of length 2r + 1 and height h = 2m + 1 . 

An example of the schematic representation of a path ojg G ^2r-2s with a given number of 
c and d marked return steps and a separating vertex V2s is shown in figure IT^ As there are no 
steps below y = 1 we have pushed the hovering path down by unit height and consider it as a 
Dyck path. 

separating vertex 




Figure 19: An example of a schematic representation of a Dyck path with a given number of c 
and d marked returns and a marked separating vertex V2s- 

As in the proof of lemma 01 inserting an up step in each path at the position of V2s and 
taking the union over all possible positions of V2s starting at the last c return and up to but 
not including the first d return (or from the beginning if there are no c returns and to the end 
if there are no d returns) replaces the pair of Bubbles on either side of the separating vertex by 
the set Ballot paths of height one and one step longer, see figure OUk-bl. Replacing each marked 
return step with an up step (and hence increasing the height of the path by 2 each time) then 
produces a Ballot path of length 2r + 1 and height 2m + 1 as required, see figure OUb ). 

□ 



Corollary 8. 



Z2r{l\l;Ki,K2) 



-{l-cd)CT 



2\2 



A'^{z'-z') 



((l + c)2-cA2)((l + d)2-dA2) 
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■* 2r + 1 *■ 

Figure 20: a) Taking the union over positions of the marked separating vertex produces a 
schematic height one Ballot path as in h), replacing each marked c and d return by an up 
step produces a Ballot path of length 2r + 1 and height 2m + 1 as in c). 

Note: This converts to the integral formula of Derrida et al ^ (BIO) with = e*^ and and using 
a contour integral to pick out the constant term. It is related to the u) expansion (see later). 
(BIO) was obtained by finding the eigenvectors of C2 and is therefore an evaluation of Z^'^\ 



Proof. With w = m. and using (|ITT^ 

Z2r-(1|1; Ki, K2) = CT 



(1 — cw){\ — dw) 

Now symmetrize the denominator by multiplying numerator and denominator by (1 — cui){l — 
dw). 

The result follows using 

(1 — uP'){l — cw){l — din) = 1 — cd — vu'^ + cduP' + {c + d){w — in) 

Because the rest of the expression is now symmetric the contribution from the last term vanishes 
by replacing w hy w and both w'^ and w'^ can be replaced by {w'^ + w'^)/2. □ 
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5 The "cj" expansion and phase diagram of the ASEP model. 



With = c/(l + c)2 corollary m may be written in the form 



Z2ril\l;Ki,K2) 

which may be expanded to give 

Z2r(l|l;«;i,K2) 



c — d 



CT 



(l-u;,A2)(l-u;rfA2) 



(5.1) 



c — d 



where 



Z2rito) = CT 



A2^'(l-z2) 



(5.2) 



(5.3) 



1 -a;A2 

The asymptotic form of Z2ri^^) as r ^ oo was obtained in & and will now be used to study 
the phase diagram for the ASEP model. First we outline the method by which the asymptotic 
form was obtained. 

Notice that expanding the factor (1 — luA'^)^^ in 1)5. 3() in powers of to gives an infinite series 
which is only valid for c < 1 which is the point at which u; as a function of c passes through its 
maximum value j. Instead we use 1)5. 3() to obtain a recurrence relation. Thus noting that 



UJ 



^2r 



l-a;A2 



A 



2r-2 



and substituting in ()5.3|) gives 



UjZ2r{(jj) = —Cr-1 + Z2r-2{'jj)- 



(5.4) 



where Cr is the Catalan number, H3.4() . Solving H5.4|) with Z2 = (1 + c)^ gives 



r-l 



Z2r{^c) = ''(1 + C - ^ C. 

j=Q 



(5.5) 



which on substituting for Uc in terms of c must give a polynomial in c. 
For the ASEP model ujc = a(l — a) and 

4uJc l-|l-2a| 



j=0 



2Wf. 



2a(l - q) 
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so 



^ oo oo 



j=0 



3=r 



I -2a) e{l-2a) + ^Cji 



where 6{.) is the unit step function. This is [H] equation (3.61). 
The asymptotic form for r ^ oo was obtained using 



vrr 



3/2 



and approximating the sum by an integral with the result 



(5.6) 
(5.7) 



l-2a 

, ,r+l 

2 4'' 

4r 



a < 2 



« > 2 



and hence 



/<(«) = i-l? 

2 4'' 



/>(«) = 4a;, 



(l0g(4^))"' 



" < 2 



" > 2 



(5.8) 



(5.9) 



As a — > |, log(4^) ~ ^^'^^a'* and these results agree with [3] equations (48)-(50). 
In the phase diagram there are therefore 

• three special regions i?i = {a > ^, /? > R2 = {a > {3, {3 < Rs = {a < (3,a < ^} 
The partition function is obtained from that in R2 by interchanging a and /?. 

• three special lines Li = {a = (3 < L2 = {P = > L3 = {a = ^, /5 > ^} The 
partition function L3 is obtained from that in L2 by interchanging a and /3. 

• a special point where the lines meet P = {a = |,/3=|} 
Table^shows the asymptotic form of Z2r-(1|1; ^i, ^2) for the above cases 
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Ri 


i?3 


P 








i?2 



a 



Figure 21: The various regions of the ASEP phase diagram. 







p . />(")-/>{/?) 




P : 4'' 











Table 1: Asymptotic forms of Z2r-(1|1; ^i, ^2) in the regions defined in the phase diagram in 
figure El 

6 Recurrence relations for the partition function and correlation 
functions. 



6.1 Recurrence relations for the partition function. 

The various formulae for G{z) when substituted in H4.16() yield recurrence relations for Zt{y\l; ki, K2) 
For example, using the identity 



1 



(1 - uJcA^){l - a;,A2) " ^ ' (1 - uj,A^){l - u^A^) ' 



1 + 
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substitution in (|5.1|) and using the CT formula (|3.4jl for Catalan numbers leads to, for r = 0, 1, 



Z2r(l|l;Ki,K2) = — T^C'r+l + (t^c + ^d)^2r+2(l|l;«:i,K2) - Wc^rf Z2r+4 (1 1 1; , K2) . (6.1) 

c — a 

Substitution in terms of ni and K2 yields, for r = 2, 3, . . . 

(k2 - l)^2r(l|l;«^l,'t2) = (k2(ki + K2) - 2Ki)Z2r-2(l| 1; , K2) + K\Z2r-A{l\l] Kl, K2) - K2Cr-l 

(6.2) 

which may be initialised by Zo(l|l; ki, K2) = 1 and Z2{1\1; ki, K2) = ki + K2- 
The following identity 

1 _^^{a + P)z/A-aPz^/A'^ 



(1 - az/A)(l - pz/A) (1 - az/ A){1 - pz/A) 

when substituted in 1)4. 28|) gives, using the CT formula H3.2() for Ballot numbers, for y = 1,2, . . . 
and odd t + y > 3 

Zt{y\l;Ki,K2) = Bt-i^y-2 + {a + (3)Zt-i{y + 1|1;ki,K2) - a(3Zt-2{y + 2\l; ki, K2) (6.3) 

which relates partition functions on lines of constant t + y. The partition function for y = 1 
is determined by (|6.1() and the following proposition then the determines Zt{2\l; ki , H2) which 
provides the initial condition for ()6.3|) . 

Proposition 3. For r = 1,2, . . . 

aZ2r-l{2\l;Ki,K2) = Z2r(l|l;Ki,K2) " ^2r (1| 1; 0, /3) 

Proof. From corollary IHl for r = 1, 2, . . . 

r m r 

Z2r{l\l; Ki , K2) = ^ i?2r-m-l,m-l ^ a*^™" *+ ^ B2r-m-l,m-lP"^ 
m=l i=l m=l 

r—1 m 

= aY, B2r-n.^2,n. J] a^^"^"' + ^2r(l|l; 0, P) 
m=0 i=0 

and the result follows from corollary with t = 2r — l,y = 2. □ 



Finally, substituting the identity 



1 + + (6.4) 



{1 - CZ'^){1 - dz-^) (l-Cz2)(l-(iz2) 
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in (|4.33jl and using (|4.16|) leads to the recurrence, for t = 0, 1, 2, ... , and odd t + y > 1 

Zt{y\l;Ki,K2) = Bt+i,y + {c + d)Zt{y + 2|1;ki,K2) - cdZt(y + 4|1; ki, Kg). (6.5) 

which relates partition functions along lines of constant t and may be initialised using the above 
relations to find the partition functions for y = 1, 2, 3 and 4. 

6.2 Recurrence relations for the correlation functions of the ASEP model. 

The probability of finding particles at positions ii,i2, ■ ■ ■ ,in is, using 

where the un- normalised n-point correlation function Gn{ii,i2, ■ ■ ■ ,in',r) is given by 

Gn{ii,i2,...,in;r) =< W\&^-^D&^-''-^D . . .&'^-'"~^"^DC''"'-\V > . (6.7) 

where C = DE. This expression may be thought of as replacing C by D in < VF|C"^'|y > at each 
of the positions ik, k = 1, . . . ,n which is equivalent to replacing the Ej matrix in a Cj = DjEj 
product by a unit matrix. Thus in the path representations Gn{ii,'i'2, ■ ■ ■ ,in',''') is obtained by 
modifying the allowed step definition such that for k = 1, . . . ,n, the step = e2i^ (beginning 
at j; = 2ik — 1 and ending at x = 2ifc) is always an up step and has weight 1. We will say that 
Sfc is a forced up step. This is illustrated in figure [ 





r V \up^ 








f \ ^'^vy V \ 


up X \y 









2ii-l 



2^2-1 2^3-1 2^4-1 



Figure 22: The n-point correlation path interpretation (n = A is shown above): each step starting 
at X = 2ifc — 1, /c = 1, . . . , n must he an up step with weight one. 



6.2.1 The case a = (3 = I 

In the case a = /? = 1, or ki = K2 = 1, it is shown in that 

Z2.(l|l; 1, 1) = < W\C'-\V > \a=f3=i = Cr+u (6.^ 
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a Catalan number. This may also be seen in terms of the third path representation (correspond- 
ing to the third matrix representation of section [21). < VF3|(L>3£'3)^|V3 > |a=/3=i is just the 
total number of One Up paths of length 2r and these biject to Dyck paths of length 2r + 2 by 
adding an up step at the beginning and a down step at the end of each path (the dotted steps 
in figure 122)1 . The result follows since number of Dyck paths of length 2p is well known to be 
the Catalan number Cp. The result also follows from 1)4. 29() by setting c = d = and noting 
that i?2r+i.i = Cr+i which is the m = term. 
It is also shown in |4|, equation (88) that 



Gnihjk, ...,in;r)= ^ Cp„Gn-l{k,i2, ■ ■ ■ 

Pn=0 



I) 



(6.9) 



This may be derived combinatorially as follows. Again we use the third path representation and 
to avoid special cases we imagine that the paths are extended to y = by a further down step. 
For each path which contributes to Gn{ii,'i2, ■ ■ ■ ^'in]f) we determine a subpath uon which starts 
with the last forced up step, s„, and ends when the path returns to the same height for the first 
time. This subpath contains a Bubble which we suppose has length 2p„ so that the subpath has 
length 2pn + 2- see figure 05k . Immediate return corresponds to pn = ^ and the maximum value 
of Pn is determined by the condition, (2i„ — 1) + (2p„ + 2) = 2r + 1, that there are no further 
steps beyond ojn- 

a) b) 

I — delete — | 





2r - 2pn - 2 



Figure 23: a) Schematic path picture for the correlation functions for a = (3 = 1, showing the 
bubble following the last forced up step, b) Removing the bubble leaves a shorter Dyck path 

We can now define a new path obtained by deleting u;„ and joining the two (possibly empty) 
resulting sub-paths which remain. This path contributes to Gn-i{ii,i2, ■ ■ ■ ,in-i',r — Pn — !)• 
The result follows by partitioning the paths contributing to Gn according to the value of For 
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a given value of p„ the number of configurations is therefore the product G„_i and the number 
of configurations of which is equal to the number, Cp„, of Dyck paths of length 2pn- 
For k = 1,2, . . . n, let 

n 

qk = n-k + ^Pj (6.10) 

j=k 

then as pointed out in [IJ, equation ()6.9() may be iterated or, combinatorially, n Bubbles may 
be removed, to give the explicit formula 

Gnihjh, ■ ■ ■ ,in]1^) = • • • Cp-^Cp2 ■ ■ ■ Cp^^Cr-q^ (6.11) 

pi>0 p„>0 

where the upper limits are Pk = r — ik — Qk+i with Qn+i = 0. 

This formula was previously conjectured by Derrida and Evans '3' on the basis of computer 
calculations up to r = 10. 

6.2.2 General a and /3. 

Equation (45) of [3] is the case n = 2 of the following proposition which we now prove using a 
lattice path representation. An algebraic proof was given in 0]. 

Proposition 4. For l<«i<i2<'<in<'' — 1 

r-i„-l 

Gn{k,i2, ■ ■ ■ ,in;r) = ^ CpGn-i{ii,i2, ■ ■ ■ ,in-i;r - p - I) 

■r-i„+l 

,l2, ■ ■ ■ , In— I ) *n 1) B2r-2i„-p,p-2P^. (6.12) 

p=2 

Notes: 

• When a = (5 = 1 this reduces to H6.9() since using 

fc+i 

B2k-p,p-2 = Ck 

p=2 

the second sum of 1)6.12(1 is just the missing term p = r — in oi the first sum. 

• In constructing a proof it was found that the first representation in terms of Jump-Paths 
was simpler to use than the third which we used in the special case of the previous section. 
It was explained after definitional, that Jump-Step paths (representing Z2r) never intersect 
y = 0. However this is not the case in calculating Gnihjk, ■ ■ ■ ^in]f) since the paths are 
modified by the forced up steps. This allows y = to be visited and then a forced up step 
returns the path to y = 1. Notice that the down step leading to y = has a weight (3. 
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Proof: Partition the modified Jump-Step configurations according as the last forced up 
step, Sn, which starts at height y = 2k, k > 1 (case A) or y = (case B) ~ see figure OH 

a) , , , , 





— delete — 








y \ 










2in 


- 1 



b) 




Figure 24: a) Case A: Last forced up step starts at y = 2k, k > 1. b) Case B: Last forced up 
step starts at y = 0. Note, the "bubbles" represent Jump Step subpaths. 



Case A: A subpath a;„ may be identified in exactly the same way as in the case « = /? = ! 
but in this case the Bubble of length 2p obtained by deleting the first and last steps of tj„ is not 
an elevated Dyck path since it may contain jump steps. However the number of configurations 
of iLin is still equal to Cp (see below) for any value of k and since u;„ avoids y = 1 its steps are 
unweighted. Thus on partitioning the paths according to the length {2p+2,p = 0, 1, . . . r — i„ — 1) 
of iOn a factor Cp may be removed from the sum over paths having the same value of p. When 
iOn is deleted the remaining steps form a weighted path of length 2r — 2p — 2 which has only 
n — 1 forced up steps. Summing over configurations of this path for given k and then summing 
over k > 1 gives Gn-iiii,i2, ■ ■ ■ , in-i', r — p — 1). The first term of the proposition formula is 
thus derived provided that we can show that the number of configurations of oJn is Cp. 

Now the paths u;„ of length 2p + 2 biject to P2p\ by vertical translation through distance 2k 
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and from (jmH) and (im 

l-^2p;ll = ^2p\a=0,i3=l = ■^2p'* Ui=0,K2=1 = (6.13) 

The last equality follows since when ki = the One Up paths which visit y = have zero weight 
and the remaining paths, which have weight 1 when K2 = Ij biject to Dyck paths by vertical 
translation through unit distance. This result is also proved combinatorially in fltJj . 

Case B: The weighted sum over paths may be factorised into three parts. 

(i) A factor which arises from the subpath consisting of the first 2i„ — 2 steps. The subpath 
ends at y = 1 and has only n — 1 forced up steps, therefore the sum over these subpaths 
yields G„_i(ii,i2, • • . ,in-i;i„ - 1). 

(ii) A factor /5 which arises from the next two steps which visit y = and return to y = 1. 
The second of these is the last forced up step Sn having weight 1. 

(iii) A factor arising from the subpath consisting of the remaining 2r — 2i„ steps which is a 
jump step path beginning and ending at y = 1, avoiding y = 0. The weighted sum over 
these paths is obtained by setting a = in the normalising factor for paths of length 
2r — 2in, thus using corollary El 

r—i 

Z2r-2i„\a=0 = -B2r-2i„-m-l, m-l/?™'- 

m=l 

The product of these three factors yields the second term of the formula. 

7 Conclusion 

We have shown that the normalisation of the ASEP can be interpreted as various lattice path 
problems. The lattice path problems can then be solved using the constant term method (CTM). 
The combinatorial nature of the CTM enables us to interpret the coefficients of the normalisa- 
tion polynomials as various un-weighted lattice path problems - usually as Ballot paths. One 
particular form has a natural interpretation as an equilibrium polymer chain adsorption model. 
The "w" form of the normalisation is particularly suited to finding the asymptotic expansion of 
normalisation and hence the phase diagram. We also formulate a combinatorial interpretation 
of the correlation functions. 
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The lattice path interpretations enable us to make connections with many other models. In 
particular, because of the strong combinatorial nature of the CTM we are able to find a new 
"canonical" lattice path representation. In a further paper |IT ] we show that this representation 
leads to an understanding of a non-equilibrium model (the ASEP model) in terms of a related 
equilibrium polymer model. Also, having extended the polymer chain model so that the end- 
points have arbitrary displacements from the surface the method of Gessel and Viennot |13j . 
jl4j (see also ) may be used to express the partition function of a network of non-intersecting 
paths as a determinant. In particular, the case of two paths gives the partition function for a 
vesicle model with a two parameter interaction with a surface. A bijection between these vesicles 
and compact percolation clusters then enables an analysis of the properties of the clusters 
attached to damp wall to be made. These applications will also be the subject of a subsequent 
publication. 

It is also of combinatorial interest to understand how the various path problems might be 
related. Clearly they are related algebraically as they are all just representations of the same 
algebra (and thus related by different similarity transformations). In a subsequent paper, |inj . 
we show combinatorially (using bijections and involutions) how the various path representations 
are combinatorially equivalent. 

8 Acknowledgements 

Financial support from the Australian Research Council is gratefully acknowledged. JWE is 
also grateful for the kind hospitality provided by the University of Melbourne and RB is also 
grateful for the kind hospitality provided by Royal Holloway College, London University. We 
would like to thank Andrew Rechnitzer and Aleks Owczarek for reading the manuscript and 
providing useful comments. 



38 



APPENDIX 



A Specimen partition functions. 



Zo(l|l;Ki,At2) =1 (A.l) 

Z2{1\1;Ki,K2) =Kl + K2 (A. 2) 

Z4(l|l; Ki, K2) ='^1^ + K2 + 2kiK2 + H2^ (A. 3) 

Zq{1\1; Ki , H2) =Kl^ + 2 K2 + 2 Ki K2 + 3 Kl^ K2 + 2 + 3 Kl K2^ + ^2^ (A.4) 
Z8(l|l; Ki, K2) =^1*^ + 5K2 + 4K1 K2 + 3ki^ K,2 

+ 4 Ki^ K2 + 5 K2^ + 6 Ki /t2^ + 6 Ki^ K2^ + 3 /t2^ + 4 Ki ^2^ + ^2^ (A.5) 

Zo(3|l;Ki,At2) =0 (A.6) 

^2(3|1;ki,aj2) =1 (A.7) 

Z4(3|1;ki,K2) =2 + Ki + K2 (A.8) 

Z6(3|l; Ki, K2) =5 + 2ki + + 3k2 + 2ki K2 + «^2^ (A. 9) 
Z8(3| 1; Ki , K2) =14 + 5 Ki + 2 + + 9 K2 + 6 Ki K2 + 3 K2 + 4 ^2^ 

+ 3kiK2'^ + K2^ (A.IO) 

Zt(l|3;Ki,K2) =/t2Zt(3|l;Ki,K2) (A.ll) 

Zo(3|3;ki,K2) =1 (A.12) 

Z2(3|3;ki,K2) =2 (A.13) 

Z4(3|3;ki,K2) =5 + K2 (A.14) 

Zq{3\3;ki,K2) =14 + 4k2 + kiK2 + K2^ (A.15) 

Z8(3|3; Ki , K2) =42 + 14 K2 + 4 ki K2 + K2 + 5 K2^ + 2 ki K2^ + K2^ (A.16) 



B Constructive proof of proposition ^ 

The formula 

Zt{y\y'; Ki,^,2) = CT[A'zy{zy' + 
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clearly satisfies the general equation (|4.7|) and the initial condition (|4.4|) provided that 

cT[zy+y^u{z)] = {). 

In order to satisfy the boundary condition (|4.11l) 



(B.l) 



CT[A*z(zJ'' + U{z)zy')] = KiCT[A'-^zizy^ + C/(z)z^')] + K2CT[A'-h\zy^ + ^/(z)/^')] 

or 

Cr[A*-2F'-i(A2 -Ki- K2zA)] = -Cr[A*-V(z)zf'-i] 

where 

V{z) = (1 - {Ri + R2)z'^ - R2z'^)U{z). 
Replacing z by z everywhere under the CT operation leaves the value unchanged so 

V{z) = -(A^ - Ki - KzK.2) = z'^{k2 -1 + {Ri + R2)z'^ - z"") 

and ()B.1|) is satisfied provided that y + > 2 or y* > 1 since y > 1 . Hence in this case y* > 1 

U{z) = z^{-l + K2G{z)). 

When y* = 1 the boundary condition is satisfied if 

Cr[A*-2(A2 -Ki- K2zk)] = CT[A*-2(z2 _2 + z'^-Ki- K2{1 + z^))] = -CT[N'-'^V {z)] 
and replacing z^ by in the second expression gives 

V{z) = R1 + R2 + R2Z^ - z^ 

which satisfies ()B.1|) and hence 

U{z) = z^{-l + G{z)). 
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